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Abstract 

We demonstrate that the Hagedorn-like growth of the number of observed 
meson states can be used to constrain the degrees of freedom of the underlying 
effective QCD string. We find that the temperature relevant for such string 
theories is not given by the usual Hagedorn value Th ~ 160 MeV, but is 
considerably higher. This resolves an apparent conflict with the results from 
a static quark-potential analysis, and suggests that conformal invariance and 
modular invariance are indeed reflected in the hadronic spectrum. We also find 
that the D± = 2 scalar string is in excellent agreement with data. 
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It is by now a well-established fact that many aspects of hadronic physics can 
be successfully modelled by strings. Indeed, an effective "QCD string" theory would 
simultaneously explain the existence of Regge trajectories, linear confinement, the 
exponential rise in hadron-state densities, and s- and t-channel duality. In fact, it 
has recently been proposed that string modular invariance might even explain relative 
meson/baryon abundances [[J. 

Over the years many different string models have been proposed for describing 
the QCD color flux tube deemed responsible for quark confinement in mesons: early 
examples include the scalar (Nambu) string, the Ramond string, and the Neveu- 
Schwarz (NS) string, and more recent examples |2j include Polyakov's "rigid string" , 
Green's "Dirichlet string", and the Polchinski-Strominger effective string. While 
all of these models are endowed with a certain number D± of bosonic degrees of 
freedom on the two-dimensional string worldsheet (corresponding to the vibrational 
and rotational string degrees of freedom in the external spacetime), the early models 
place additional, purely internal degrees of freedom on the string worldsheet. The 
later string models instead introduce effective interactions for the spacetime degrees 
of freedom which alter their short-distance behavior. 

Conformal invariance nevertheless plays an important role in each case. In the 
early models, conformal invariance is exact (though classical), and indeed an infinite 
array of models of this type can be constructed by placing on the string worldsheet 
virtually any two-dimensional conformal field theory (CFT). The quantum excitations 
of the corresponding new worldsheet fields should then also produce new hadronic 
states. Conformal invariance in some of the more recent models, by contrast, is an 
effective symmetry, valid in the long-string limit. The effective CFT's which emerge 
in this limit are nevertheless equally responsible for reproducing the salient features 
of hadronic spectroscopy. 

A comparison with data is therefore necessary in order to constrain the class of 
strings appropriate for modelling hadronic physics, or, more generally, to test the va- 
lidity of the string approach by determining the extent to which conformal invariance 
is actually reflected in the observed hadron spectrum. In particular, we shall use data 
to constrain the worldsheet central charge c of the effective QCD string, since this 
parameter is a unique measure of the degrees of freedom in a general two-dimensional 
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theory. Most previous efforts in this direction have employed the static quark poten- 
tial method, and favor values c ~ 2. We propose using the well-known Hagedorn-like 
growth in the number of meson states as an independent method constraining not 
only the total central charge, but also its distribution between spacetime and inter- 
nal degrees of freedom. While the value Th ~ 160 MeV commonly taken for the 
Hagedorn temperature is too low to agree with the above results (yielding c ~ 7), 
we will show that for comparisons with strings an alternative treatment is necessary. 
This will resolve the conflict between these two methods, and in so doing provide ex- 
perimental evidence that the QCD string possesses conformal - and indeed modular 
- invariance. As a by-product, we will also find that the D± = 2 scalar string is in 
excellent agreement with data. 

Let us first review the basic ideas behind the static potential V(R) between two 
quarks a distance R apart. Modelling confinement as a color flux tube, one can show 
[|J that this potential takes the exact form 



where a is the string tension of the flux tube and M 2 < is a constant independent 
of R. While the first (linear) term in the large-i? expansion of V(R) represents 
the classical energy in the effective string, the second "pseudo-Coulomb" term is 
an attractive universal quantum correction (or Casimir energy) which arises due to 
transverse zero-point vibrations of the string. As such, this term is to be distinguished 
from the true attractive Coulomb term which has the same form and which arises 
at short distances from gluon exchange. The form of the exact string result ([!]) 
indicates that while aR plays the role of a string "momentum," the quantity Mq 
appears as the string "rest mass" (or ground-state energy) j5j . The fact that M 2 < 
(or equivalently that the long-distance pseudo- Coulomb term is attractive) implies 
that the ground state of the effective QCD string is tachyonic, yet this causes no 
inconsistency in the large- R limit 0. 

If we assume the dynamics of the color flux tube to be modelled by a two- 
dimensional conformal field theory, then this ground state energy M and the central 



V(R) 



[aRf + M, 2 ]' 
aR + (M 2 /2a) iT 1 + £>(iT 2 ) , 



(1) 
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charge c of the corresponding worldsheet theory are related by 

a' M 2 = h-c/24 (2) 

where a' = (2ira)~ 1 is the Regge slope, and h is the conformal dimension of that 
primary field in the worldsheet theory which produces the ground state. In most 
cases this primary field is merely the identity field with h = 0, so that the coefficient 
of the pseudo- Coulomb term directly yields the corresponding central charge. In all 
other cases, however, the coefficient of this term yields information concerning only 
the difference h — c/24 = — c/24. This is dramatically illustrated in the Ramond 
string: here h — c/24 — (D — 2)/16, whereupon c = 0, the ground state is massless, 
and the long-range pseudo-Coulomb term is absent. 

By fitting the parameters of Eq. (|T|) to heavy-quark spectroscopic data and/or 
results from lattice QCD, many authors have attempted to determine the ground- 
state energy Mq 2 and thereby the central charge of the QCD string. While the string 
tension a is generally found to be in good agreement with the Regge-trajectory value 
a! ~ 0.85 (GeV) -2 , values of c have been obtained throughout the range < c < 4, 
clustering near c ~ 2. Perhaps the largest source of error in these methods is the fact 
that they rely upon a full separation of the effects of the long-range Coulomb term 
from those of the true Coulomb interaction: while the latter are a priori unrelated 
to c, fits to spectroscopic data and/or lattice QCD results undoubtedly contain their 
contributions. Furthermore, as discussed, c is not always a true measure of the 
degrees of freedom in the string worldsheet theory. 

There exists, however, a different approach towards determining the central charge 
c, one which avoids all of the above difficulties and is complementary to that involving 
the static quark potential. In string theory (or more generally in any CFT), the 
number or degeneracy of states g n at any excitation level n is given by the coefficients 
in a certain polynomial Xh(%) called the character of the sector [h] of the worldsheet 
CFT: 

oo 

Xh (x) = x-^Y,9nX n . (3) 

ra=0 

Since the spacetime mass M n of a given string excitation level n is given by a'M n 2 = 
n — c/24, we see that the ground state energy M 2 in each sector [h] is indeed given by 
Eq. (Q). Thus, the static quark potential method, by fitting M , is essentially a test 
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of the string n — > limit. However, as has been well-known from the earliest days of 
the dual-resonance models, information can also be extracted from the high-energy 
limit (n — > oo), for in this limit string theories predict an exponential rise in the 
degeneracy of states g n with excitation number n: 

g n ~ A [C 2 (n - 5/24)] ~ B e c ^ n -~ c/2i 

= A (M/T H )~ 2B e M ' TH . (4) 

Here A, B, and C are constants, with Th = (CVa')~ 1 . The form of these expressions 
demonstrates that Th is the famous Hagedorn temperature 13], a critical tempera- 
ture beyond which the partition functions of such theories (and indeed all of their 
thermodynamic quantities) cannot be defined. The exponent B also has profound 
physical consequences. Since the internal energy of such a hadronic system near Th 
grows as U{T) ~ (Th — T) 2B ~ 3 [|J, we see that Tjj is a true maximum temperature 
if B < 3/2, and merely the site of a second-order QCD phase transition otherwise. 

What makes these observations useful for our purposes, however, is the fact that 
Th in Eq. (|J) is directly related to the total central charge of the underlying CFT: 

Indeed, this result holds independently of the sector [h] in which our string states are 
presumed to reside, yielding a value for the true central charge c rather than c and 
providing a test independent of the static quark potential method. Furthermore, the 
exponent B can also be interpreted in terms of an underlying string theory, for B is 
universally related to D±, the effective number of spacetime dimensions for transverse 
string oscillations (or equivalently the number of uncompactified bosonic degrees of 
freedom in the worldsheet CFT ||): 

B = \{3 + D ± ) . (6) 

Thus, a fit of Eq. (§) to hadronic data yields information concerning not only the total 
number of degrees of freedom, but also their effective distribution between spacetime 
and internal excitations. 

In practice, however, a number of subtleties arise which must be addressed before 
an adequate fit can be performed. First, by tabulating the experimentally measured 
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masses Mj and widths r, of observed mesons JED], we have calculated the density 
p cxp (M) of meson states, where p exp (M) = Y,iW i S{M;M i ,Y i ). Here WJ = 7«(2ii + 
l)(2Jj + l) is the number of states per resonance (with j$ — 1 for charge self-conjugate 
states and 7^ = 2 otherwise), and 5 represents a statistical distribution function. 
The result is plotted in Fig. 1, with error bars determined by varying S between 
Breit-Wigner, Gaussian, and fixed-width (Tj = 200 MeV) distributions. These errors 
also include an estimate of the uncertainties resulting from the possible subtraction 
of quark masses. While this density clearly experiences the predicted exponential 
growth over much of the plotted mass range, we see that the rate of growth sharply 
diminishes beyond 1.7 GeV. This is attributable to experimental difficulties, for at 
higher energies it becomes harder to distinguish mesons from background. Indeed, a 
comparison with explicit quark- model calculations ]lT] shows 1.7 GeV to be the first 
energy where fewer mesons are observed than predicted. At the other extreme, the 
data below 0.3 GeV depends on pion contributions whose small masses result from 
chiral symmetries which string theory is not expected to model. Thus, we shall limit 
our attention to the experimental meson data for 0.3 < M < 1.7 GeV. 

This in turn requires a more sophisticated treatment of the theoretical string 
predictions than was sketched above. In particular, taking a' ~ 0.85 (GeV) -2 , we 
see that a mass near 1.7 GeV corresponds only to a string excitation of n « 3. While 
the n — > 00 asymptotic function quoted in Eq. (^) is remarkably accurate even for 
relatively small values of n, it differs from the true string degeneracies g n for n < 3 
by as much as 95% for the scalar string. We thus require from string theory a more 
precise functional form, one which includes a sufficient number of subleading terms 
so that the true values of g n for the known strings are accurately reproduced. It is 
a straightforward matter [E| to determine the first set of these subleading terms, 
however, and together these yield the following improved "asymptotic" form for the 
state degeneracies gM and corresponding density p(M): 



g M ~ V2tt A CI\u\(0 where i = M/T H 

= (2a , M)- 1 P stnng(M) . (7) 

Here L u \ is the modified Bessel function of order \u\, with v = 1/2 — 2B; note that 
this result reproduces Eq. (||) in the limit M — > 00 jn|. By substituting the proper 



values of A, B, and Th for known strings, we have verified that Eq. (|7j) is indeed 
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accurate over the required range to within 2%. This is fortunate, since the forms of 
any additional subleading terms are dependent on model-specific parameters other 
than B and Tjj. 

We then performed a fit comparing the experimental and theoretical values of 
J™. p(M)dM as a function of m, with = 0.3 GeV. Our results are as follows. 
Taking B = 5/4 (i.e., = 2) resulted in a best fit with 

Tjj = 300 MeV =>■ c = 1.97 , (8) 

where we have taken a' = 0.85 (GeV) -2 . This is remarkably close to the central 
charge c = 2 of the D± = 2 scalar string, and demonstrates that a string picture 
is indeed consistent with the data obtained from counting the numbers of hadronic 
states. Note, in this regard, that the original 1967 fits by Hagedorn yielded 
Tjj « 160 MeV, or c ~ 7, in clear contradiction with the results from the static 
inter-quark potential. However, the fundamental difference is the functional form 
to which the fits are made, for Hagedorn's bootstrap-motivated functional form || 
p(M) = a(M 2 + /i 2 )~ 5//4 e M//TH with ii 500 MeV has no connection to the string- 
theory result in Eq. fl7|) and implicitly assumes B = 7/4 in the M — > oo limit. Indeed, 
it is only upon taking the string-motivated functional form in Eq. (|7|) that agreement 
is obtained. 

In light of interesting recent proposals ]14| for non-scalar QCD strings with values 
B 7^ 5/4, it is nevertheless important to place general constraints in (B,Tjj) param- 
eter space. In Fig. 2, the singly- and doubly-shaded regions indicate the (B,Tjj) 
values for which fits to the data with a xVd.o.f. < 1 can be obtained: the dashed 
line indicates the central value of the fits (i.e., the "best fits"), and the heavy solid 
lines indicate the border of the region allowed by string theory (corresponding to 



the constraints B > 3/4, c > D± fL5f ). We see that there indeed exists a region of 
overlap between the data-allowed and string-allowed regions; furthermore, the "best 
fit" line passes directly through this overlap region in the range 3/4<5<5/4. In 
the case of scalar strings for general D± (which lie along the upper curved boundary 
of this string region), we see that the "best fit" line intersects this scalar string line 
almost exactly at the expected value B = 5/4, or D± = 2. 

The overlap region can be further narrowed if we constrain the normalization 
constant A in Eq. (|7|). Within the allowed region in Fig. 2, we find that A^ t ranges 
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from 0(1) to > O(10 3 ). String theory, however, places the precise limit v4 str ing < 
y/2ir (4:71o:'Th 2 ) u , with equality occurring for the scalar string. Thus, if we assume 
that there are 36 independent strings which contribute to the (it, d, s)-quark meson 
spectrum (corresponding to 36 quark degrees of freedom: 9 possible quark/anti- 
quark flavor combinations, 4 spin states, and one color singlet state), we can require 
r = Afit/A^^g < 36. This then restricts us to the lower (doubly-shaded) portion of 
the data-allowed region in Fig. 2. In Fig. 3 we have plotted versus (c, D±) that region 
of allowed parameter space which satisfies all three of these string constraints. The 
best fit line is superimposed. 

We have already seen in Eq. (|8|) that the c = D±_ = 2 scalar string lies directly on 
the best fit line; we now see from Fig. 2 that this point is also exactly on the r = 36 
border. This implies that the scalar string should accurately model the absolute 
numbers of states in the meson spectrum, and not merely their rate of growth. In 
Fig. 1, for example, we have superimposed on the meson data the actual numbers 
of states predicted by 36 copies of the scalar string (solid line). The agreement is 
excellent. 

In conclusion, we have seen that the density of meson states is consistent with 
string-theoretic predictions; moreover, estimations of the central charge of the QCD 
string obtained via measurements of the appropriately-defined Hagedorn tempera- 
ture are now consistent with those independently obtained via static quark potential 
methods. This latter agreement is especially significant, for these two methods de- 
pend separately on quantities which are a priori independent: the rate of growth 
of the numbers of mesons, and the energy of the ground state of the corresponding 
flux tube. Indeed, as originally noticed in Ref. only an underlying effective 
two-dimensional conformal invariance — in particular, modular invariance and the 
associated symmetry under x — > — 1/x in Eq. @ — serve to relate them. Our results 
thus constitute strong additional evidence that the confinement phase of QCD is 
consistent with an effective string theory in which conformal symmetry and modular 
invariance play a significant role. This certainly warrants further study. 

Acknowledgments 

We are pleased to thank A. Mironov for initial discussions, and M. Li and R. 



8 



Myers for comments on the manuscript. This work was supported in part by NSERC 
(Canada) and FCAR (Quebec). 

References 

[I] P.G.O. Freund and J.L. Rosner, Phys. Rev. Lett. 68, 765 (1992); J.R. Cudell 
and K.R. Dienes, Phys. Rev. Lett. 69, 1324 (1992). 

[2] A.M. Polyakov, Nucl. Phys. B268, 406 (1986); M.B. Green, Phys. Lett. B266, 
325 (1991); J. Polchinski and A. Strominger, Phys. Rev. Lett. 67, 1681 (1991). 

[3] R. Hagedorn, Nuovo Cimento 56A, 1027 (1968). 

[4] O. Alvarez, Phys. Rev. D24, 440 (1981); M. Ltischer, K. Symanzik, and P. Weisz, 
Nucl. Phys. B173, 365 (1980). 

[5] L. Brink and H.B. Nielsen, Phys. Lett. 45B, 332 (1973). 

[6] P. Olesen, Phys. Lett. 160B, 144 (1985). 

[7] See, e.g., Ref. || and references cited therein; also A. A. Bykov, A.V. Leonidov, 
and A.D. Mironov, Mod. Phys. Lett. A4, 125 (1989). 

[8] S. Frautschi, Phys. Rev. D3, 2821 (1971); K. Huang and S. Weinberg, Phys. Rev. 
Lett. 25, 895 (1970); N. Cabibbo and G. Parisi, Phys. Lett. 59B, 67 (1975). 

[9] The general relation is 5 = 3/4 — k/2 where k is the modular weight of the 
characters x associated with the worldsheet CFT ||12|| . Each transverse spacetime 
dimension (and the Liouville mode, if present) contributes k = —1/2. There 
are, however, no contributions from compactified bosonic worldsheet fields (or 
equivalently from the worldsheet fermions of the R/NS strings). 

[10] Particle Data Group, K. Hikasa et a/., Phys. Rev. D45, SI (1992). 

[II] S. Godfrey and N. Isgur, Phys. Rev. D32, 189 (1985). 



9 



[12] G.H. Hardy and S. Ramanujan, Proc. Lon. Math. Soc. 17, 75 (1918); I. Kani 
and C. Vafa, Commun. Math. Phys. 130, 529 (1990). Since I- V {x) = I u (x) if 
(/ G Z or a; > 1, we restrict ourselves to Bessel functions with positive order 
in Eq. (0); these are better-behaved in the x ~ 0(1) and unphysical (v ^ Z) 
regions. 

[13] Since the states in a non-interacting string theory are infinitely narrow and 
populate discrete energy levels M n , their density is p(M) = qm J2n S(M — M n ) = 
gM&ct'M) J2 n 3(. a 'M 2 + c/24 — n). We replace the latter sum of ^-functions with 
unity in order to provide a more realistic estimate of the generic broadening 
effects that interactions would have on the string spectrum. Furthermore, in the 
asymptotic limit n — * oo, this yields p(M) = guj 'AM where AM is the mass 
difference between adjacent string levels. 

[14] M. Caselle, R. Fiore, F. Gliozzi, P. Provero, and S. Vinti, Phys. Lett. B272, 272 
(1991). 

[15] It is possible, however, to violate the constraint c > D± through renormalization 
group effects: see, e.g., A.B. Zamolodchikov, JETP Lett. 43, 730 (1986). 

[16] L. Brink and H.B. Nielsen, Phys. Lett. 43B, 319 (1973); W. Nahm, Nucl. Phys. 
B114, 174 (1976). 



10 



Fig. 1: Number of meson states with masses < M as function of M, compared with 
scalar-string result. 
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Fig. 2: Data- and string-allowed values of (B,Th), as discussed in text. Best data 
fits (dashed line) and D± = 2 scalar string point (dot) also shown. 
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Fig. 3: Values of (D±, c) satisfying both data and string constraints, with best fit 
line superimposed and scalar string shown (dot). 
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